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The photonic band dispersion and density of states (DOS) are calculated for the three-dimensional 
(3D) hexagonal structure corresponding to a distributed Bragg reflector patterned with a 2D tri- 
angular lattice of circular holes. Results for the Si/Si02 and GaAs/AlGaAs systems determine the 
optimal parameters for which a gap in the 2D plane occurs and overlaps the ID gap of the multilayer. 
The DOS is considerably reduced in correspondence with the overlap of 2D and ID gaps. Also, the 
local density of states (i.e., the DOS weighted with the squared electric field at a given point) has 
strong variations depending on the position. Both results imply substantial changes of spontaneous 
emission rates and patterns for a local emitter embedded in the structure and make this system 
attractive for the fabrication of a 3D photonic crystal with controlled radiative properties. 
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I. INTRODUCTION 

Photonic crystals are being intensively studied with 
the goal of achieving a photonic band gap as well as con- 
trol of light propagation and light-matter interaction at 
visible frequencies A complete photonic gap in 

three dimensions (3D) has been demonstrated for the di- 
amond lattice of dielectric spheres J 5[ , the fee lattice of 
air spheres (or "inverse opal") |||7j]7 the "yablonovite" 
H , the "woodpile" |J, and other more complex 3D ge- 
ometries ftLOUlll . All these structures are difficult to fab- 



ricate at optical wavelengths: they often require the use 
of a bottom-up procedure, like e.g. in the case of in- 
verse opals, where the template is first built by self- 
assembling of dielectric spheres in a solution and the 
voids are subsequently filled with a high refractive in- 
dex material |^-|l4|| . As an alternative, two-dimensional 
(2D) photonic crystals can be fabricated with a top-down 
approach based on lithography and etching. This pro- 
cedure can be realized at sub-micrometric wavelength 
scales and allows for the controlled introduction of lin- 
ear and point defects. A main problem however concerns 
the control of light propagation in the third (vertical) 
dimension: even in the case of a waveguide-embedded 
2D photonic crystal, radiation losses due to out-of-plane 
diffraction in the vertical direction cannot be eliminated 
for photonic modes which lie above the light line [jl5] |l7| . 

In this work we explore theoretically another possi- 
bility for achieving control of light propagation in 3D, 
namely the use of distributed Bragg reflectors (DBRs) 
patterned with a 2D lattice in the layer planes. DBRs 
are dielectric multilayers which are periodic in one di- 
mension (ID) 18-20), i.e., they represent ID photonic 
crystals and have a band gap for propagation of light 
along the multilayer axis. By patterning and etching a 



DBR with a 2D lattice which possesses a photonic gap, a 
3D photonic crystal with uniaxial (or biaxial) symmetry 
is obtained: if the ID gap of the multilayer is made to co- 
incide with the gap of the 2D structure, a common band 
gap for light propagating along the main crystal axes can 
be achieved. Although a complete photonic band gap in 
3D is not expected (due to the angular dependence of the 
ID and 2D gaps, or else to the anisotropy of the Brillouin 
zone), an appreciable control of the photonic density of 
states (DOS) and consequently of the spontaneous emis- 
sion properties may be realized. In such a situation, the 
etched DBR represents a favorable system for the study 
of spontaneous emission control in 3D photonic systems, 
since local (and possibly active) defects can be introduced 
at the level of ID epitaxial growth and of 2D lithographic 
design. 

Here we focus on the simplest and most promising 
structure, namely a DBR patterned with a triangular 
lattice of holes: the planar structure is known to dis- 
play a 2D photonic gap common to both polarizations 
of light, if the dielectric constant and air filling fractions 
are large enough p|,pl|-p3]| . We calculate the 3D photonic 
band dispersion and the photonic DOS for parameters 
which are representative of the Si/Si02 system (a typical 
high-index contrast DBR) and of the GaAs/Al^Gai-ajAs 
system (with low index contrast). In particular, we de- 
termine the conditions for the ID and 2D photonic gaps 
to overlap: since the ID stop band of the DBR changes 
upon patterning, and the 2D lattice of holes is formed 
in a stratified medium, the overlap of the two band gaps 
represents a sort of selfconsistency problem. The density 
of states is found to be considerably reduced in corre- 
spondence with the overlapping 2D and ID gaps. Also, 
the local density of states (which is the relevant quan- 
tity for spontaneous emission rates) shows strong varia- 
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tions and additional reductions as a function of position. 
These findings show that sizeable changes of spontaneous 
emission rates and patterns of a local emitter embed- 
ded in this anisotropic 3D structure should occur. The 
present results may also serve as guidelines for experi- 
mental groups which are attempting the fabrication of 
these structures. 



II. THEORETICAL METHOD 




The geometry of the assumed structure is shown in 
Fig. la. The z axis is taken along the DBR growth di- 
rection. We denote by h,h the thicknesses of the two 
DBR layers and by E\ , £2 their respective dielectric con- 
stants. The 2D pattern in the xy plane is taken to be 
a triangular lattice of pitch a, with a basis in the unit 
cell consisting of circular rods of radius r and dielectric 
constant Eh- In the present paper we consider only air 
holes (eh = 1). The 3D Bravais lattice is simple hexago- 
nal and the corresponding Brillouin zone with the main 
symmetry points is shown in Fig. lb. The DBR period 
L = li + l 2 is taken as a fixed length scale of the prob- 
lem, in the sense that the photonic frequencies are given 
in dimensionless units ujL/(2ttc): indeed, if such a struc- 
ture is grown, a single epitaxial DBR can be patterned 
with various lateral periods at fixed L. The 3D photonic 
lattice is fully determined by the ratios h/L, a/L, r/a as 
well as by the dielectric constants £1, £2, £/i- 

In order to calculate the photonic band structure we 
start from the second-order vector equation for the mag- 
netic field, 



Vx— VxH=— H, (1) 
e(r) c l 

and use the well-known plane-wave expansion method. 
When N reciprocal lattice vectors G are kept in the ex- 
pansion, the equation is transformed into a 2N x 2N ma- 
trix eigenvalue problem which contains the Fourier trans- 
form £q q/ °f the inverse dielectric constant e ~ 1 (r ) . This 

is evaluated by the method of Ho, Chan and Soukoulis 
||, which consists of calculating the Fourier transform 
£q q' of e(r) and inverting numerically the resulting 
matrix. This procedure yields a rapid numerical conver- 
gence when an energy cutoff is imposed for the number 
of reciprocal lattice vectors. 

The density of states is defined as 

P H^^^- Wn (k)), (2) 

n k 



(b) 




FIG. 1. (a) Schematic structure of a DBR patterned with a 
triangular lattice of circular holes; (b) the hexagonal Brillouin 
zone with the irreducible wedge. The origin of the coordinates 
in (a) is taken at the center of a hole in the middle of layer 1. 

where n is the index of photonic bands and the sum 
over wavevectors extends over the first Brillouin zone. 
The local DOS is similarly defined as 

P (w,t) = E E l E «k(r)| 2 ^ " w„(k)). (3) 

n k 

In evaluating Eq. || the integral of the squared electric 
field is normalized to the unit cell volume, therefore DOS 
and local DOS have the same dimensions and can be di- 
rectly compared. In order to calculate both quantities we 
adopt the linear tetrahedron method p^-p6[. The DOS 
is already converged with a mesh of about 400 k-points in 
the irreducible wedge of the Brillouin zone (see Fig. lb). 
Calculating the local DOS requires a mesh which spans 
the whole Brillouin zone. 

We now calculate the Fourier transform of the dielec- 
tric tensor, defined by 




for the patterned DBR structure (V is the unit cell vol- 
ume). 
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FIG. 2. Bands (a) and density of states (b) for a A/4 Si/Si02 multilayer. Parameters: ei = 12, £2 = 2, h/L = 0.3, a/L = 1, 
and r/a — 0.45. The dashed line in (b) represents the photon DOS in an average isotropic medium, while the dotted line is the 
DOS of a uniaxial medium, whose dielectric tensor components are deduced from the photonic bands in the long-wavelength 
limit (see text). 



We write the 3D reciprocal lattice vectors as G = 
(Gu,G z ), where G|| is the projection in the xy plane; 
similarly, r = (r||,z). The unit cell volume is written as 
V = AL, where A is the unit cell area for the 2D lattice. 
For a general patterning, the Fourier transform can be 
expressed in terms of a 2D structure factor 



A ./did 



iGii -r 11 



dr 



(5) 



where the integral extends only over the unpatterned (or 
dielectric) region. Note that F(Gy = 0) = / is the 2D 
dielectric filling fraction. For G = 0, the 3D Fourier 
transform is simply the average dielectric constant 



,h 



e(G|, = 0, G z = 0) = (-^ + -±£ 2 )/ + e h (l - f) = e 



L 



L 



(6) 



For G 7^ 0, it can be evaluated by subtracting the vanish- 
ing quantity y- J coll e lGr dr from Eq. (Q): the integrand 
e(r) — £h now vanishes in the holes and the integral over 
the dielectric region yields the 2D structure factor (||), 
with the result 



e(G|| ^0,G Z =0) = (| £l 



£2 -£ h )F(G||), (7) 



Note that the last expression is valid for any G|| (=0 or 
7^ 0). The above equations (@)-(@) hold for a patterning 



with a general 2D lattice and holes of arbitrary shape: 
the 2D structure factor can be written as 



nGii)=« Gl .o/- (1-^.0)1/^ 



e iG ii- r n rfr||. (9) 



For the present case of a triangular lattice of circular 
holes, the unit cell area is A = \Z3a 2 /2 and the 2D struc- 
ture factor becomes 



F(G||=0) = 1-— , (10) 



F(G|, 0) = J^r), (11) 

where J\ is the Bessel function of first order. 



III. RESULTS 

A. Si/Si0 2 system 

We first consider the case of Si/SiC>2 multilayers, mod- 
eled by the dielectric constants £\ = 12 and £2 = 2. The 
Si/SiC>2 system has the appealing features that the di- 
electric contrast is large, implying that a high rejection 
ratio in the stop band is obtained with a small number 
of periods, and that the etching of the two materials is 
a technologically well controlled process. Another high- 
index contrast system described by similar parameters is 
GaAs/Alox (oxidized AlAs) §7|. 
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In Fig. 2 we show the photonic bands and density of 
states in the case of a Si layer width l-y/L = 0.3: this is 
close to the A/4 condition for the multilayer (i.e., when 
the two layer thicknesses are inversely proportional to 
their refractive indices) and it maximizes the ID gap 
along the vertical direction. For the 2D lattice we choose 
a/L = 1 and a hole radius r/a — 0.45: the latter value 
produces a full photonic band gap in the pure 2D case, 
when the refractive index contrast between dielectric and 
air is strong enough j|j2l|-^]. Inspection of Fig. 2a shows 
that a ID gap along the T-A direction is indeed realized 
around ojL/(2itc) — 0.4, but also that no 2D gap oc- 
curs in the 2D projection of the Brillouin zone (r-M-K-r 
directions). The reason is that the average dielectric con- 
stant of the multilayer is too low to support a 2D band 
gap. In order to achieve a 2D photonic gap, the A/4 con- 
dition must be abandoned and multilayers with a higher 
Si fraction should be considered. 

Concerning the photonic density of states, Fig. 2b 
shows that the DOS increases like uj 2 at small frequencies 
(like for free photons) but pronounced structures occur 
when the photonic band dispersion deviates from linear- 
ity. In this and the following figures for the DOS, the 
dashed line represents the photon DOS in a homogeneous 
and isotropic medium with the same average dielectric 
constant e av of the etched DBR (Eq. ^), while the dotted 
line is the DOS of a homogeneous but uniaxial medium, 
whose dielectric tensor components are obtained from the 
in-plane and T-A dispersions at small wavevector p8[ . 
The dotted curve agrees with the numerically calculated 
DOS in the low-frequency (long-wavelength) limit; the 
dashed curve represents a reference DOS of an "average" 
medium, and is useful in order to tell whether the DOS 
of the etched DBR is reduced or increased with respect 
to that of the average medium. 

In Fig. 3 we show the photonic bands and DOS of 
three Si/Si02 multilayers with l\/L = 0.8 and a/L = 
0.8, 1.2, 1.8: in all these cases the hole radius r/a = 0.45. 
A full 2D band gap in the r-M-K-r plane can be recog- 
nized in Fig. 3c (around ujL/(2ttc) = 0.4) and in Fig. 3e 
(around ljL/(2ttc) = 0.27), since the average dielectric 
constant of the multilayer is now close to the Si value. In 
Fig. 3a, the 2D gap should coincide with the second-order 
ID gap at the F point: however, the 2D gap is actually 
closed by a photonic mode which starts at the lower edge 
of the ID gap at about ujL/(2ttc) = 0.51. On increasing 
the ratio a/L the 2D gap opens and decreases in energy, 
and for a/L — 1.8 (Fig. 3e) it overlaps the first-order ID 
gap at the A point. 

In the case shown in Fig. 3e, it would seem that the ID 
gap along T-A is closed by a photonic mode which starts 
at the lower 2D gap edge around uL/(2irc) — 0.27. How- 
ever, since this mode is nondegenerate, it does not have 
the same symmetry of the electromagnetic field (which 
belongs to a twofold degenerate representation for any 



wavevector k along T-A) and it cannot couple to an ex- 
ternal beam incident along the T-A direction: this mode 
is "symmetry uncoupled" |4|,|2^-|3l| and it cannot have 
any observable effect on the reflectance for a plane wave 
incident along the multilayer axis. Indeed, calculations of 
the optical properties [3^ j indicate that the etched DBR 
still behaves as a ID reflector with a well-defined stop 
band along T-A and with a reflectivity which tends to 
unity as the number of periods increases. 

Concerning the photonic DOS, most features are simi- 
lar to those of Fig. 2b already commented: however it is 
interesting to observe that a pronounced minimum occurs 
in correspondence with the 2D gap along r-M-K-r (e.g., 
around ujL/{2itc) — 0.27 in Fig. 3f). In the frequency 
window of the 2D gap the DOS is reduced not only with 
respect to the "average" DOS (dashed line), but also in 
comparison with the DOS of the long wavelength limit 
(dotted line). The reduction of the DOS is particularly 
pronounced in the cases of Figs. 3 (b) and (f), when the 
2D gap overlaps the ID gap. Therefore, a reduction of 
spontaneous emission rates and a strong change of emis- 
sion pattern is expected under these conditions: for the 
parameters of Fig. 3 (e,f), since spontaneous emission 
is suppressed in the xy plane and along the vertical di- 
rection (due again to symmetry mismatch between the 
emitted photon and the nondegenerate mode in the gap) , 
the emission of a point-like dipole will be preferentially 
directed along the diagonals of the 3D structure. 

In Fig. [| we show the local density of states (LDOS) 
for the "optimal" parameters of Fig. 3(e,f), i.e., when 
the 2D and first-order ID gap overlap. The LDOS is 
evaluated at three different positions along a line con- 
necting two neighboring holes in the midplane of layer 1 . 
The general behavior of the LDOS is similar that of the 
DOS, with a u> 2 increase at small frequencies, a reduction 
in correspondence with the common 2D- ID gap and pro- 
nounced structures at higher frequencies. As discussed in 
connection with Eqs. |2|,|[, the DOS and LDOS are chosen 
to have the same units and can be directly compared, 
or in other words, the spatial average of the LDOS over 
all positions gives again the DOS. It can be seen that 
the LDOS is increased over the DOS of Fig. 3(f) for a 
position inside a hole (Fig. [|a or b), while it is consider- 
ably reduced at the center of the narrow dielectric region 
between two holes (Fig. |]c). The combined effects of a 
common 2D-1D gaps and of a change in position (i.e., of 
the electric field profile) results in a sizeable reduction of 
the LDOS: at r = (0,0.5a,0) and ujL/{2ttc) = 0.28, the 
LDOS is equal to 4.4, compared to about 33 for the av- 
erage medium. Therefore, the numerical results confirm 
that a strong reduction of spontaneous emission rates can 
be achieved for an emitter embedded in an etched DBR, 
when the parameters of the structure are properly cho- 
sen. 
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FIG. 3. Bands (left) and density of states (right) for Si/SiCh multilayers. Parameters: ei = 12, £2 = 2, h/L = 0.8, and 
r/L = 0.45. Figs. (a,b): a/L = 0.8, Figs. (c,d) a/L = 1.2, Figs. (e,f) a/L = 1.8. Dashed and dotted lines in the right panels: 
same as in FigM 
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FIG. 4. Local density of states (LDOS) for Si/Si02 multilayers. Parameters: e 1 — 12, £2 = 2, h/L = 0.8, a/L = 1.8, 
and r/a — 0.45. (a): LDOS at position (0,0,0), i.e., at the center of a hole in the middle of layer 1; (b): LDOS at position 
(0, 0.3a, 0); (c): LDOS at position (0, 0.5a, 0), i.e., at the center of the dielectric vein between two neighboring holes. 



B. GaAs/AlGaAs system 



IV. CONCLUSIONS 



We now consider the case of GaAs/AlGaAs multilay- 
ers, modeled by the parameters e\ = 11.7 and £2 = 10.5: 
the two values correspond to the dielectric constants of 
GaAs and AlGaAs, respectively, at 1 eV p3| ]. Since the 
two refractive indices are large and close to each other, 
the photonic energies in dimensionless units do not de- 
pend sensitively on the layer thicknesses and we 
choose them to be equal, lx = h = 0.5 L: this is also 
close to the A/4 condition for the multilayer. In Fig. ^ 
we show the photonic bands and DOS with these param- 
eters, assuming two different values a/L = 1.2 and 1.8 
and taking again a large hole radius r/a = 0.45. A 2D 
gap along T-M-K-r is formed (around luL/(2ttc) = 0.36 
in Fig. |p, around ujL/(2ttc) — 0.24 in Fig. ||c) and in 
the case of a/L = 1.8 it overlaps the small first-order 
ID gap at the A point. Like in the case of Fig. 3e, the 
nondegenerate mode starting at the lower edge of the 2D 
gap is not coupled to a plane wave incident along the 
T-A direction and it will not affect the normal-incidence 
reflectivity. Concerning the density of states, there is 
again an increase like u> 2 at small frequencies followed 
by several structures. A reduction of the DOS occurs in 
correspondence with the 2D gap, with a minimum which 
is somewhat less pronounced than for the etched Si/Si02 
multilayers: this is due to the smaller refractive index 
contrast between GaAs and AlGaAs. Changes of spon- 
taneous emission rates and patterns will also occur in the 
GaAs/ AlGaAs system, but to a lesser extent than for a 
high index contrast structure. 



Systematic calculations of photonic bands, total and 
local DOS have been undertaken for the hexagonal pho- 
tonic structure corresponding to a distributed Bragg re- 
flector patterned with a triangular lattice of holes. A few 
representative examples have been shown for the case of 
Si/Si02 and GaAs/ AlGaAs multilayers, focusing on the 
most favorable situations for achieving a full 2D photonic 
gap in the T-M-K-r plane. The parameters for which 
the 2D gap overlaps the ID gap formed along the multi- 
layer axis have been determined. A nondegenerate mode 
whose dispersion along T-A overlaps the ID gap is sym- 
metry uncoupled from a plane wave propagating along 
T-A and it will not affect the normal-incidence reflectiv- 
ity 

The conditions for the overlap of a 2D and a ID gap 
are rather restrictive and require a careful structure de- 
sign: the average refractive index of the multilayer has 
to be large (thus the layer thicknesses must not obey 
the A/4 condition) and the air fraction of the 2D lattice 
must be such that a full 2D gap develops. The overlap 
of the 2D with the first-order ID gap is found to occur 
around a/L ~ 1.8 for the investigated structure, which 
is favourable in view of patterning the multilayers with 
lithography and etching techniques. 

The photonic density of states has a pronounced min- 
imum in correspondence with the common 2D-1D gap, 
particularly for the case of the Si/Si02 system. The local 
DOS depends strongly on the position and is considerably 
reduced at the center of the dielectric veins between two 
holes. The combined effects of overlapping 2D- ID gaps 
and of the position dependent DOS will yield a strong re- 
duction of spontaneous emission rates and considerable 
changes of emission patterns, with a redistribution of the 
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emission mainly along the diagonal direction of the 3D 
structure. 

Specific calculations of the optical properties of etched 
DBRs (reflectivity and spontaneous emission) are under- 
way pi] . The anisotropic structure made of the etched 
DBR is concluded to be an interesting possibility for the 
realization of 3D photonic crystals where a controlled in- 
troduction of defects is possible and whose spontaneous 
emission properties can be tailored to a large extent. 
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FIG. 5. Bands (left) and density of states (right) for GaAs/AlGaAs multilayers. Parameters: ei = 11.7, £2 = 10.5, h/L = 0.5, 
r/a = 0.45. Figs. (a,b): a/L = 1.2, Figs. (c,d) a/L = 1.8. Dashed and dotted lines in the right panels: same as in Fig.B. 
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